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I demonstrate that 4 He adsorbates in carbon nanotube materials can be treated as one- 
dimensional interacting gas of spinless bosons for temperatures below 8 K and for coverages such 
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adsorbate-adsorbate interactions are studied within the scheme of virial expansion approach. The 
theoretical predictions for the specific heat of the interacting adsorbed gas are given. 
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I. INTRODUCTION 

The adsorption of gases in nanotube-based materials has been recently a subject of considerable interest and many 
theoretical and experimental studies focused on this phenomenon have been reported. 1-12 The interest in the subject 
stems partially from the possibility to use these materials as efficient gas containers for hydrogen storage. 13 Another 
cause of the interest is that the nanotube materials provide a very specific potential energy environment for the gas 
atoms and molecules. In particular, the nature of this environment is such that it reduces the effective dimensionality 
of adsorbates 8,9 which at sufficiently low temperatures behave as a one-dimensional (ID) gas. This provides an 
excellent opportunity to study the interactions in the ID gas. The problem of N particles interacting mutually via 
binary interaction potentials in one dimension has been thoroughly investigated in the literature and there exist 
exact quantum and classical solutions for very specific functional forms of the interaction potential. 14-17 The aim 
of this article is to investigate a realistic system in which the adsorbate atoms (molecules) interact with a relatively 
complicated binary potential that is attractive at large and repulsive at short interadsorbate separations. 18 

The outline of this article is as follows. In Sec. II, the behavior of single atom adsorbed on the surface of a 
bundle of single-wall carbon nanotubes (SWCNT) is discussed. The range of temperatures in which the isolated 
adsorbates exhibit effective one-dimensional behavior is discussed for 4 He atoms. In Sec. Ill the calculation of 
second virial coefficient for the interacting gas in ID is briefly outlined. A fully quantal approach is followed as the 
gas of interest is composed of 4 He atoms for which quantum effects are essential. In Sec. IV the specific heat of 
adsorbed He gas is predicted and the effects of He- He interactions are discussed. The results obtained are compared 
with those for exactly solvable models 15,16 and qualitative agreement is found. The calculation of specific heat in 
combination with experiments, some of them quite recently reported 19,20 , is expected to yield an additional insight 
in the thermodynamics of the adsorbed gas. In Sec. V, the influence of the corrugation of the nanotube (which is 
neglected in Sees. II, III and IV) and sample inhomogeneities on the thermodynamics of adsorbed gas is thoroughly 
discussed. Sec. VI concludes the article and summarizes the main results of the present work. 

II. BEHAVIOR OF ISOLATED 4 He ATOMS ADSORBED ON THE SURFACE OF A BUNDLE OF SWCNT's 

Although there is still an ongoing discussion 1-4,9,10 concerning the preferable adsorption sites for He atoms in 
SWCNT's materials, the experimental information 1-3 combined with the theoretical considerations 9-11 suggests that 
individual He atoms are predominantly adsorbed in the groove positions on the SWCNT's bundle surface. If the 
number of He atoms adsorbed in the sample is very large, then one can expect that the He atoms will also occupy 
other positions on the bundle surface. This point will be discussed in Sec. III. 

Quantum states of 4 He atoms adsorbed in the groove positions of an infinitely long bundle made of (10,10) SWCNT's 
have been discussed in Ref. 9. It was found that the low-energy part of the 4 He excitation spectrum exhibits a typical 
ID behavior with characteristic \j\f~E singularities present in the density of states, g{E). The density of states 
does not exhibit gaps which is a consequence of the neglect of the corrugation of the carbon nanotube. In this 
approximation there are no potential barriers for the adsorbate motion along the groove. The severeness of this 
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approximation and its influence on the results to be presented shall be discussed in Sec. V. In Fig. 1, the low-energy 
part of the 4 Hc density of states per unit length of the groove is presented, which was calculated as described in Ref. 
9, i.e. the three-dimensional Schrodingcr equation was numerically solved to yield the complete set of bound states. 
As can be seen in Fig. 1, all the excitations with energies between -22.7 meV and -18.88 meV pertain to essentially 
one-dimensional 4 He atoms. The excitations in this regime of energies represent the activated, free motion of 4 Hc 
atoms along the groove characterized by a ID wave vector K y (y-axis is oriented along the groove), as discussed in 
Ref. 9. The transverse profile (in xz plane, perpendicular to the groove direction) of the 4 He wave function is the 
same for all these excitations and can be represented by a narrow, Gaussian-like 2D function (see panel A of Fig. 4 
in Ref. 9). 

At -18.88 meV, another band of states becomes available to the isolated adsorbates. In this band, the transverse 
profile of the 4 He wave function is different from the ground-state profile (see panel B of Fig. 4 in Ref. 9). As discussed 
in Ref. 9, the population of higher (excited) bands causes transition from the effectively ID behavior of 4 He atoms to 
2D, and eventually 3D behavior. 

For the purposes of this work, it is sufficient to note that the separation between the lowest ID band and the first 
excited band is quite large (3.82 meV) which immediately suggest that the higher bands are poorly populated in a 
significant range of temperatures. The width of this temperature range can be evaluated from the known density of 
states. The total density of states can be represented as a sum of the lowest band density of states, go(E) and the 
density of states representing all other transverse excitations, g a , 

g(E)=g (E)+g a (E). (1) 

The lowest band density of states per unit length of the groove is given by 9 

where the mass of the adsorbate is to, = —22.7 meV is the ground state energy, and 6 is the Heaviside function. 
The total number of adsorbates, N, is given by 

/>oo 

N = L g(E)f(E,T)dE, (3) 
where L is the total length of the groove, and the Bose-Einstein distribution function is given by 

f(E,T)= * . (4) 

Here, ks is the Boltzmann constant, T is the temperature and v is the chemical potential. The number of adsorbates 
in the lowest band, No, can be calculated as 

/•OO 

N = L g (E)f(E,T)dE, (5) 

once the chemical potential has been determined from Eq. (3). 

In Fig. 2, I plot the ratio N /N as a function of temperature and for three different linear densities of adsorbates, 
n, n = N/L. From this figure one can conclude that the noninteracting 4 He gas can be treated as effectively being ID 
for temperatures smaller than about 8 K (13 K) since for these temperatures more than 99 % (95 %) of 4 He atoms 
occupy the lowest energy band. 

III. VIRIAL EXPANSION APPROACH TO TREAT THE INTERACTING GAS IN ID 



The approach to be presented here assumes that all the adsorbed atoms are in the groove positions on the bundle 
surface. Thus, the approach can not be applied to the situations where the number of adsorbates is so large that other 
positions on the bundle surface become occupied by the adsorbates. One could envisage the situation where a very 
large number of atoms is adsorbed on the bundle surface. From this phase, one could get to the phase where all the 
adsorbates are exclusively in the groove positions by desorbing all the atoms which are not in groove positions. These 
atoms are more weakly bound than those in the grooves 9 and will desorb from the sample at lower temperatures. This 
fact enables the experimental realization of the ID phase of interest to this work. Similar arguments can be applied 
to ID heavy adsorbate phases studied in Ref. 12. 
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The virial expansion approach to treat the imperfect (interacting) quantum gas is well known 21-23 and in this 
section it will be only briefly outlined with reference to a quantum gas in ID. 

The basic idea of the virial expansion approach is to represent the so-called gas spreading pressure, <f), as a power 
series of the gas density, n, 

oo 

t34> = Y J Bi{L,p)n l . (6) 
1=1 

Here, /3 = l/fceT and coefficients Bi are the virial coefficients. The virial coefficients can be obtained by comparing 
the expansion of the gas spreading pressure in the power series of fugacity, z = exp(pv), 

(34>=-\nQ{z,p,L) = Y J h{L,p)z l (7) 
i=i 

with the expansion in Eq. (6). 6/ is the l-th cluster integral obtained as the coefficient in the power series expansion 
of the logarithm of the grand partition function, Q, in terms of fugacity. Since the grand partition function is given 
by 

OO 

Q{z,L,(3) = Y,Zi{L,(3)z l , (8) 

1=0 

where Zi are the quantum partition functions for I particles, the formulas for cluster integrals can be obtained by 
comparing the expansions in Eqs. (7) and (8). The general form of bi is given in Ref. 21. For the first two cluster 
integrals one has 

i z x 
h = T 

h = ^£2. (9) 

Eliminating z from equations (6) and (7) by expressing it in terms of linear density ra, yields the relation between 
the virial coefficients and cluster integrals. Using Eq. (9), one can obtain the relations between virial coefficients and 
quantum partition functions. Explicitly, for the first and second virial coefficient one has 



B l = 1 



1 z. 



B * = L {-2-z!)> ^ 
The quantum expression for the partition function of N particles in ID is 

Zn = dy 1 ...dy N ^2^>* a (y 1 , ...,y n )exp[-PH( Pl , ...,p N , yi , ...,y N )}^ a ( yi , ...,y n ), (11) 

J a 

where pi and yi, i — 1,...,N represent ID momenta and coordinates of the N gas particles, respectively. The 
dynamics of N gas particles is described by the hamiltonian, H(pi, ...,pn, yi, Vn)- A complete set of quantum 
states describing N particles is denoted by {a}. The wave functions, ty a arc assumed to be properly normalized 
and symmetrized according to the statistics satisfied by the gas particles. It is easy to show 21 " 23 that the partition 
function for one ID spinlcss particle is given by 

Zi = j, (12) 



where A is the thermal wavelength, A = y2irh 2 /3/m. Since I am going to consider 4 He atoms, the consideration of 
spinlcss particles will suffice. As shown by the authors of Ref. 21, the spin degrees of freedom can be considered 
(if needed) after the spinless problem has been solved. In the problem of spinlcss particles (anti)symmetrization of 
the wave function is performed solely in the coordinate space. All expressions which follow do not consider the spin 
degrees of freedom. 
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For the calculation of the second virial coefficient, B 2 , given by Eq. (10), one needs to calculate the partition 
function for two interacting particles. This is an easy task when the corrugation of the nanotubes can be neglected, 
since in that case the interacting two-body problem can be reduced to free motion of the center of mass and the 
relative motion representing a particle of reduced mass /1 = m/2 in an external potential 21- . The motion of the 
center of mass can be represented by a wave function for the free particle of mass 2m, exp(ik cm Y) /\/L, where k cm is 
the wave vector of the center of mass motion and Y is the center of mass coordinate, Y — (yi + 1/2) The relative 
motion can be described by a wave function of relative coordinate, y = y\ — y 2 , which is denoted by £, c (y), and which 
satisfies ID Schrodinger equation, 

~^ + v(\y\)\uy) = ^Uy)- (is) 

The binary potential representing an interaction between the two gas particles is denoted by v(\y\) = v(\y\ — 2/2 1 ) - The 
set of "relative" quantum states is denoted by {c}. This set consists of a finite number of bound states denoted by 
{&}, and the continuum of states which can be numbered according to the wave vector k associated with the motion 
of the particle with reduced mass in the region where the interaction vanishes (large y). The energy of quantum state 
|c) is denoted by e c . Relative wave functions behave in the asymptotic regime (y — > 00) as 

Uy) -> (14) 

The ID phase shifts, which are nonvanishing due to the presence of the interaction potential v, are denoted by f](k). 

The partition function for the two noninteracting particles [w(|y|) = 0] in ID can be calculated without invoking Y 
and y coordinates. Its form results solely from the requirement of the (anti)symmetrization of the total wave function. 
Explicitly, 

z " -- ± UT\> (15) 

where the uppcr(lowcr) sign is for spinless boson(fcrmion) gas. The second virial coefficient for the noninteracting gas 
is thus 

< = ^- < 16 > 

The superscripts (0) in Eqs. (15) and (16) indicate that the expressions correspond to the noninteracting quantum 
gas. As often noted in the literature 22,23 , the finite value of coefficient for noninteracting quantum gas reflects 
the so-called statistical attraction for bosons (negative ) and statistical repulsion for fermions (positive B^ , see 
Eq. (6)). The expression for B^ should be compared with Eq. (3.15) of Rcf. 21 which pertains to noninteracting 
gas in two dimensions. 

The second virial coefficient can be calculated as 

B 2 = B<°> + AV2^ [exp(-/34°)) - exp(-/3e c )] , (17) 

c 

where is the set of "relative" energies for two noninteracting particles. Eq. (17) has been obtained by performing 
integration over the center of mass coordinate in the expression for partition function of the two interacting particles, 
Z 2 (Eq. (11)). 

The formula for B 2 which is convenient for numerical implementation can be obtained by replacing the summation 
over {c} in Eq. (17) with two summations, one going over the bound states, {&}, and the other over the continuum 
states, {k}. One can pass from the sum over states {k} to the integral over wave vector k by introducing the density 
of states in k space, which can be obtained from the ID phase shifts, w(k). 21 I finally obtain 

ft 03) = 5f(/3) - AV2][>xp(-/^) - ^ ^exp dk, (18) 

where the dependence of virial coefficient on the temperature (or (3) is emphasized. This formula is very similar to 
the one obtained for 2D and 3D gas, although the numerical factors (such as y/2) and units of B 2 (in ID, B 2 has 
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units of length since the linear density has units of inverse length) are different, depending on a dimensionality of the 
problem. It should also be noted that since the problem involves only one dimension, one does not obtain azimuthal 
quantum numbers which occur in the treatments of interacting gas in 2D and 3D as a consequence of the central 
symmetry of the binary potential. 

The evaluation of Eq. (18) requires the calculation of ID phase shifts which depend on the binary potential, v. 
The interaction between the two He atoms in the otherwise empty space is known to a great precision 18 . However, 
the effective interaction between the two He atoms positioned in the vicinity of a third polarizable body is different 
from the free space He- He interaction 25 . In the case of interest to this work, the two He atoms are surrounded by two 
SWCNT's and the polarization induced in the SWNCT's will modify the He-He interaction. While the polarization 
induced effects on the binary potential can be calculated for atoms physisorbed on crystalline surface 25 , the analogous 
calculation for the very specific geometry of the nanotube bundle is certainly more difficult. It is interesting to note 
here, that Vidali and Cole 26 found that the measurements of specific heat of He overlayers on graphite 27 can be more 
accurately reproduced by the effective He-He potential which is 15 % shallower from the free-space He-He potential. 
They attribute this effect to the screening of He- He interaction by the substrate. The treatment of Vidali and Cole 26 
was also based on quantum virial expansion. In another study 28 , more related to the system considered here, the 
authors found that the interaction between two He atoms adsorbed in the interstitial channels of SWCNT's has a well 
depth which is 28 % shallower with respect to the free-space interaction. The groove adsorption represents a situation 
which is "somewhere in between" the adsorption on planar graphite and in SWCNT's interstitial channels. 

In the following calculations, the He-He interaction will be described by free-space potential suggested recently by 
Janzen and Aziz 18 , but I shall also consider the scaled potential obtained from the free-space interaction by simple 
multiplication with a factor of 0.785. Thus, the scaled potential has a well depth which is 21.5 % smaller from the well 
depth of the free-space potential. This number was obtained as a simple arithmetic mean of the well depth reductions 
found for adsorption on planar graphite 26 and in interstitial channels of SWCNT's 28 . The assumed reduction of the 
well depth is quite close to the numerical estimate in Ref. 28 (24 %). The exact value of the scaling factor used should 
not be taken too seriously because the substrate induced contributions to the potential cannot be modelled by a simple 
scaling of the free-space potential. 25 The scaled potential was introduced simply to examine the effects of the details 
of interaction potential on the thermodynamics of adsorbed gas. Additionally, to obtain the effective potential in ID, 
the 3D potential should be averaged over the 2D cross-sections of the adsorbate probability density 26,29 . However, 
as the cross-section of the lowest band states is rather small, and in the light of the uncertainties of the substrate 
mediated forces, such a procedure has not been performed. 

SAPT1 potential supports one weakly bound state in ID, representing a 4 He dimer with an energy of eo = — 0.16/xeV. 
This state, being so weakly bound, is extremely extended in relative coordinate. 30,31 . The bound state energy in ID 
is significantly smaller from the one obtained by Siddon and Schick in a 2D treatment 21 , which is in accord with 
existing literature 30 . The scaled SAPT1 potential does not support bound states. The ID phase shifts and their 
derivatives with respect to relative wave vector were calculated by numerically solving the Schrodingcr equation using 
the algorithm 32 quite similar to the one described in Ref. 33. 

In Fig. 3, the calculated values of second virial coefficient are presented. The full line corresponds to calculation 
with the He-He potential suggested in Ref. 18 (SAPT1), while the dashed line corresponds to calculation using the 
scaled SAPT1 potential. The behavior of the second virial coefficient with temperature is qualitatively similar to the 
one obtained for 4 He adsorbates on graphite (2D problem) in Ref. 21. There is, however, one important difference. 
The ideal gas term, i3°, reflecting the purely quantum effect of gas statistics, decays with temperature as yT/T in ID, 
and as 1/T in 2D. 21 . Thus, the approach of second virial coefficient to its classical value as the temperature increases 
is slower in ID than in 2D. In the inset of Fig. 3 the derivatives of the phase-shifts, dr](k)/dk are plotted as a function 
of relative wave vector k. Note that the phase-shift derivatives become negative and nearly constant for large relative 
energies (wave vectors). This is a consequence of a strongly repulsive potential at short distances (hard core). Note 
also that the phase shifts of the two potentials are very different for small wave vectors, and thus one could expect that 
the two potentials produce quite different thermodynamical quantities. However, the very different behavior of the 
phase shifts is a consequence of the fact that the SAPT1 potential supports a weakly bound state whereas the scaled 
SAPT1 potential does not. Thus, in the evaluation of Eq. (18) one has to properly account for the bound state which 
exists in the case of SAPT1 potential. This "extra" term for SAPT1 potential makes the thermodynamical quantities 
derived from the two potentials quite similar, although the derivatives of the phase shifts are very different. This 
fact has been discussed for the interacting gas in 2D 34,35 in connection with Levinson's theorem 24,37 which relates 
the phase shift at zero momentum to the number of bound states. It was found 34 that a proper account of both the 
continuum and bound states eliminates discontinuities in thermodynamic properties whenever an extra bound state 
appears with the small change of the parameters of interaction potential. In the present calculation similar effect is 
found in one dimension. 

IV. SPECIFIC HEAT OF ADSORBED 4 He GAS 
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The specific heat of interacting quantum gas can be calculated from the set of virial coefficients 22,23 . I assume that 
the dominant contribution to the specific heat comes from the second virial coefficient. Thus, the results are applicable 
to a restricted range of adsorbate concentrations and temperatures. The range of a validity of this approximation can 
be estimated from a calculation of higher virial coefficients, which is a difficult task, or from direct comparison with 
experiments, as has been done for 4 He adsorbates on graphite in Refs. 21 and 36. Experiments dealing primarily with 
the specific heat of adsorbates in carbon nanotube materials have not been reported yet and those which detected 
the signature of the adsorbed gas in the overall specific heat of the sample were focused on the specific heat of clean 
nanotube materials 19,20 . 

The isosteric specific heat is given as 21 

Nk B 2 P dp 2 ' 1 j 

The second term in Eq. (19) can be calculated from Eq. (18) as 



^=W2 



^ + So + Io +/3{Si + h) _ 02{S2 + l2) 



(20) 



where 



and 



S n = ^e£exp(-/3e 6 ), (21) 

b 



k^^L^dk. (22) 



dk \ 2/i 



The expression for the term in Eq. (20) representing a deviation of the specific heat from its ideal value (where the 
interadsorbate interactions are neglected) is different from the corresponding expression one would obtain in a 2D 
treatment. The 2D expression 21 does not contain a term proportional to A and independent of the interadsorbate 
interaction potential (1/16 in Eq. (20)). The reason for this is that in 2D, the noninteracting value of B 2 , B^ is 
proportional to (3 (B^ = — A 2 /4) 21 , and therefore, its second derivative with respect to [3 vanishes. In ID case, 

oc (see Eq. 16) and d 2 B^ /d(3 2 oc /3~ 3 / 2 . This fact alone suggests that the specific heats of dilute, interacting 
boson gas in ID and 2D may be qualitatively different. 

In Fig. 4 I plot the quantity —0 2 d 2 Bij d(i 2 . The full (dashed) line represents the calculation with SAPT1 (scaled 
SAPT1) potential. Obviously, both potentials produce similar deviations. The calculation with scaled SAPT1 po- 
tential yields somewhat smaller effects of interactions on the specific heat, which is plausible since the scaled SAPT1 
potential is weaker than SAPT1 potential. An important observation is that the deviation of specific heat from its 
ideal value is negative, for ID boson gas, i.e. the inclusion of interactions reduces the specific heat of the adsorbed ID 
gas. For 2D spinless boson gas considered in Ref. 21, the deviation was found to be positive. The observed difference 
between ID and 2D results is a consequence of both the nonvanishing second derivative of the ideal term [B^\f3)] 
with respect to f3, and nonexistence of azimuthal degrees of freedom in ID treatment of the problem. 

Finally, in Fig. 5, 1 plot the specific heat of the interacting 4 He gas adsorbed in grooves of SWCNT bundles for three 
different linear densities. The full (dashed) lines represent the calculation with SAPT1 (scaled SAPT1) potential. 
Obviously, for denser gas, the specific heat is more strongly influenced by the interactions, and deviates more from the 
ideal (noninteracting) ID value (compare with Fig. 6 of Ref. 9). It can also be inferred from Fig. 5 that the description 
of thermodynamics of the adsorbed gas in terms of second virial coefficient only breaks down at very low temperatures 
where such an approach yields negative specific heat. These temperatures can be considered as the lower limits for 
the application of the present approach. For higher densities, the lower temperature limits obviously increase, in 
agreement with Eq. (19), and the presented virial expansion approach breaks down at higher temperatures. Note 
that the upper temperature for which the specific heat was calculated is 6 K. Above that temperature the higher 
bands start to contribute to the specific heat as discussed in Sec. II and Ref. 9. Although the number of particles in 
higher bands is very small at 6 K (see Fig. 1), the derivative of internal energy with respect to temperature (specific 
heat) is very sensitive even to very small occupations of the higher bands and this causes the increase of the specific 
heat observed for a very dilute gas. 

It is now of interest to compare the obtained results with those of exactly solvable ID many-particle models. In 
particular, for the gas of bosons in one dimension interacting via repulsive <5-function potential, Lieb and Liniger 
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have demonstrated that the that the energy spectrum of such a gas is identical with the spectrum of noninteracting 
Fermi gas. 15 ' 16 The same was found in the model of Girardeau 1 for ID bosons interacting with the binary hard-core 
potentials of finite radius. The Fermi energy, Ep of the corresponding ID Fermi gas is given as 

% 2 k 2 &V 2 

= ~2m^ = ~2m ' (23) 

where kp = ttN/L — im is the Fermi wave vector. This mapping (interacting Bose gas - noninteracting Fermi gas) 
allows us to easily predict the specific heat of the impenetrable ID Bose gas. Thus, at temperatures much smaller 
than the Fermi temperature, Tp = Ep/ks, 

Cy(T«T F )ocT, (24) 

while at temperatures much larger than the Fermi temperature, specific heat reduces to its classical equipartition 
value, 

CV(T>T F ) 1 

Nk~ B 2' {2b) 

The linear temperature dependence of specific heat [Eq. (24)] can be also interpreted as a signature of long- wavelength 
compression waves of the ID Bose gas (sound), and is thus also expected in a gas of ID bosons interacting with more 
complex forces. This is a collective effect and is obviously outside of the scope of the second order virial expansion, 
treating only the two-body collisions. On the other hand, the classical equipartition value of specific heat [Eq. (25)] 
is obtained also by the present approach at high temperatures [Eq. (19) and Fig. 4]. The temperature separating 
the two characteristic behaviors in Eqs. (24) and (25) is roughly given by the Fermi temperature of the equivalent 
noninteracting Fermi gas, which is proportional to the square of linear density [Eq.(23)]. Thus, for denser Bose gases, 
the interval of temperatures in which the specific heat is significantly smaller than 1/2 is larger. This is again in 
agreement with the results presented in Fig. 5. 

V. THE EFFECT OF THE CORRUGATION OF CARBON NANOTUBE AND SAMPLE INHOMOGENEITIES 



The results presented in this article are based on an idealized representation of the SWCNT materials. In particular, 
it is assumed that SWCNT bundles are infinitely long, straight and smooth which is certainly not the case in real 
materials and the bundles wiggle on large length scales. The discreteness of the nanotube results in the corrugation 
of the potential experienced by adsorbates which is not accounted for in the present approach. The influence of 
the corrugation of the potential confining the He atom to the interstitial channel of the bundle composed of (10,10) 
carbon nanotubes was studied in Rcf. 11. A model potential was used which enabled an easy examination of effects 
induced by the potential details. Substantial effect of the corrugation on the density of states of single He atom was 
predicted. In particular, for the lowest band of states, a mass enhancement factor of 2.37 was calculated. Another 
study 6 predicted a mass enhancement factor of 1.3. A large difference between these factors predicted by the two 
studies can be attributed to the larger separation of the tubes (by 0.1 A) adopted in Ref. 6 with respect to Ref. 
11. For comparison, the mass enhancement factor for 4 He adsorbed on graphite was found to be only 1.06. 26 For an 
interstitial channel surrounded by three (17,0) nanotubes, the authors of Rcf. 6 found that a change in the intertube 
separation by 0.1 A changes the effective mass of 4 He by a factor of two! 

It is quite remarkable that a small change in the intertube separation results in a very large change of the band 
structure of adsorbates. A similar sensitivity of the adsorbate bound states on the interaction potential was found in 
Ref. 9 for the groove adsorption. The details of the interaction potential are "magnified" in the spatially restricted 
regions of the interstitial channel and the groove, and the band structure of adsorbates is much more affected by the 
potential details when compared with the adsorption on planar graphite. Thus, in order to assess the effect of the 
corrugation on the calculations presented in this article, one would need to know the intertube separation distance to 
a precision better than at least 0.05 A. Note that the groove region is at the surface of the bundle where the relaxation 
effects can be expected and the separation between the tubes surrounding the groove needs not to be the same as the 
separation between the two tubes in the interior of the bundle. Additionally, one would have to know the He - single 
tube potential with a great accuracy. Even then, the calculation of the corrugation of the potential would require 
knowledge of the alignment of the two tubes surrounding the groove. It is very likely that this arrangement is not 
the same for all the grooves on the bundle surface. Furthermore, if the two tubes surrounding a groove are not of the 
same symmetry, the corrugations of the two tubes are not necessarily commensurate and the total potential for the 
groove adsorption can not be written as a Fourier series. 

All the mentioned complications are not present for the adsorption on planar graphite. In that case, a fairly reliable 
potential can be constructed 38 and written as a 2D Fouries series. However, even the inclusion of periodic corrugation 
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in the formalism of quantum virial expansion is not straightforward. The basic reason for this is that the two body 
Schrodinger equation in the textured potential background no longer separates into the two equations describing the 
motion of the center of mass and the relative motion. For the adsorption of He on graphite, Guo and Bruch 29 devised 
a perturbation treatment in which the two-body cluster integral is written as an expansion in powers of the Fourier 
amplitudes of the atom-substrate potential. They also presented the results for the second virial coefficient in which 
the corrugation effects were treated up to the second order of the perturbation series. Another study dealing with 
the problem 39 started from the tight-binding Hamiltonian and He atoms localized on particular adsorption sites. A 
series of approximations and simplifications was needed to obtain the formula for the second virial coefficient which 
is amenable to evaluation. 

It is clear from the discussion in this section that a reliable calculation of the second virial coefficient with the 
effects of the corrugation included is not possible at present, mainly because the relevant potential is not known 
with sufficient precision. Additionally, aperiodic corrugation has not been treated in the literature in this context. 
However, a simple estimate of the corrugation effects is possible if we represent the adsorbate as a particle with the 
effective mass. Such a representation of the corrugation effects is adequate for excitations with small velocity i.e. 
for the states located around the center of the ID Brillouin zone. The mass enhancement for the groove adsorption 
can be expected to be somewhere in between 1.06 (planar graphite) and 2.37 (interstitial channel). The quantity 
—(3 2 d 2 B 2 /d{3 2 was calculated using the Eqs. (20), (21) and (22), the scaled SAPT1 potential, and the effective He 
mass equal to M* = 1.3M He =5.2 amu. At T =1K (5K), its value was found to be -1.69 A(-0.52 A). This should be 
compared with the values obtained in Fig. 4. In particular, at T =1K (5K) the values obtained using the free mass of 
4 Hc are -1.92 A(-0.77 A). Thus, the reflection of the corrugation effects on the specific heat can be significant, but the 
overall trends, at least in the model of renormalized mass, are the same. For strongly corrugated potentials, it may 
be more sensible to use the band width for the characterization of the effects of corrugation rather than the curvature 
of dispersion curves at small wave vectors. 40 

The influence of other bundles on the atoms adsorbed in a groove of a particular bundle is neglected. This 
approximation obviously breaks down if two bundles touch each other. At these points, the adsorbates move in a 
potential very much influenced by both bundles in question, which may be significantly different from the potential 
of a single, infinitely long and straight bundle I used in the calculations. For the adsorption of 4 He on graphite 27 ' 36 , 
it was argued 41 that the presence of long-range inhomogencities in the graphite may act as a trap and induce the 
Bose-Einstein condensation of the adsorbed gas. The same experiments were later explained in terms of the virial 
expansion approach 21 by assuming a perfect graphite substrate. Nevertheless, a possible substantial influence of the 
sample inhomogencities on the thermodynamics of the adsorbed gas cannot be a priori ruled out. This is expected 
to be more important at low temperatures where the adsorbate thermal wavelength is large and the long range order 
of the sample may influence the adsorbate gas thermodynamics. It is hoped that the low-temperature measurements 
of the specific heat may give answers to these questions. 

VI. SUMMARY AND CONCLUSION 

4 He gas adsorbed in the grooves of single wall carbon nanotube bundles has been treated as an interacting Bosc 
gas in one dimension. It was found that this approximation should be very accurate for all temperatures below 
8 K. The interactions in the adsorbed gas are treated via the quantum virial expansion approach and the second 
virial coefficient for the interacting gas was calculated. This information was used to calculate the specific heat of 
adsorbed 4 He gas which was shown to be substantially influenced by interadsorbate interactions already at relatively 
low adsorbate linear densities (0.033 1/A). A qualitative agreement between the results obtained in this article and 
those of exactly solvable ID models is demonstrated. 
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FIG. 1. Density of states per unit length of the groove of single 4 He atom adsorbed in the groove of a bundle made of (10,10) 
SWCNT's. The low energy portion of the density of states is displayed. 



FIG. 2. Ratio of the number of 4 He atoms in the lowest ID band and the total number of 4 He atoms (No/N) as a function of 
temperature and for three different linear densities. Thick full line: N/L— 0.01 1/A. Thick dashed line: N/L=0.1 1/A. Thick 
dotted line: N/L— 0.2 1/A. Two thin dotted lines represent the 0.99 and 0.95 values of the ratio. 



FIG. 3. Second virial coefficient for ID 4 He gas as a function of temperature. Full line: Quantum calculation with SAPT1 
potential. Dashed line: Quantum calculation with scaled SAPT1 potential (see text). Inset: Phase shift derivatives, dr/(k)/dk, 
corresponding to SAPT1 potential (full line) and scaled SAPT1 potential (dashed line). 



FIG. 4. Deviation of the specific heat per unit linear density (—f3 2 d 2 B2/d[3 2 ) from its ideal value, Eq. (19). Full line: 
Calculation with SAPT1 potential. Dashed line: Calculation with scaled SAPT1 potential. 
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FIG. 5. Specific heat of 4 He gas adsorbed in the grooves of SWCNT bundles for three different linear densities (0.033 1/A, 
0.1 1/A, and 0.2 1/A). Full line: Calculation with SAPT1 potential. Dashed line: Calculation with scaled SAPT1 potential. 
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